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Abstract: The decomposition of the nucleon spin into that of quarks and gluons is related 
to twist-two GPDs according to Ji sum rule. Further decomposition into the helicity 
and the orbital angular momentum inevitably requires twist-three GPDs. In this paper 
we derive exact relations between twist-three GPDs and the canonical orbital angular 
momentum density of quarks and gluons, and check their consistency with the longitudinal 
spin sum rule. Our work demonstrates that the complete decomposition of the nucleon 
spin fits well with the framework of perturbative QCD. 
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1. Introduction 

In a by now classic paper B], Ji derived sum rules which relate certain moments of gener- 
alized parton distributions (GPD) to quarks' and gluons' individual contributions to the 
nucleon spin ~ : 

J q = 1 y dxx{H q (x) + E q (x)) , J9 = ^ J dx (H g (x) + E g (x)) . (1.1) 

Here, J 9 (J 9 ) is the fotaZ spin of quarks (gluons), namely, it is the sum of the helicity and 
the orbital angular momentum (OAM). The corresponding GPDs, H q ^ g and E q>g , are twist- 
too, and hence represent the leading contributions to cross sections in exclusive processes. 1 
In the case of J q , one can also extract the OAM of quarks L q in a longitudinally polarized 
nucleon by subtracting from J q the quark helicity ^AS which can be determined from 
independent measurements. Thus, in Ji's framework the nucleon spin decomposition reads 

- = Ji + J9 = IaS + L q + J 9 . (1.2) 
2 2 y J 

Theoretical cleanness and experimental feasibility of this decomposition scheme has boosted 
the study of GPDs to the forefront of present-day research on nucleon structure. 



However, (1.1) does not exhaust the whole content of the longitudinal spin structure. 
It turns out that one can access more detailed, interesting information if one goes to twist- 
three. An early indication of this was seen in a parton model analysis P, 0] in which a par- 



ticular twist-three GPD was shown to be related to L q . [This will be reproduced in ( 2.22; ) 



below in full QCD.] More recently, it has been recognized that a complete decomposition 
of the nucleon spin may be achieved at the twist-three level ||, |j, and this simultaneously 
resolves the longstanding controversy over whether such a decomposition is possible at all 
in QCD. [See the recent intense discussions in §, g, g g 0, |ll|, 0, [y| ||, |l|.] 



Our normalization of the gluon GPD is such that, in the forward limit, H g (x) = xG(x) where G(x) is 
the gluon distribution function. 
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Let us explicate this last point which has actually constituted the main motivation for 
the present study. Recall the well-known Jaffe-Manohar decomposition [fl6f | 

l -= l -^ + Ll an + &G + Ll an , (1.3) 

where AG is the gluon helicity, and Lean and Lean are the canonical OAMs of quarks and 
gluons, respectively. They feature the canonical momentum (ordinary derivative) x x id 
as opposed to Ji's dynamical OAM L q featuring the covariant derivative x x iD. Because 
of this, each term of the decomposition ( |l.3[) is not gauge invariant (excepting AS) once 
the interaction is turned on. Nevertheless, starting from ( |1.3|) and interpreting the fields 
as given in a particular, but arbitrary gauge, one can make the gauge invariant exten- 
sion (GIE) of it H 13 1 . The first such attempt, starting from the light-cone (LC) gauge 



expression, was made by Bashinsky and Jaffe 17] which was recently rediscovered in a 
different line of argument |10[| . Another GIE based on the Coulomb gauge was suggested 
by the work of Chen et al. (6|, [?J which also prescribed a general procedure to construct a 
GIE. Operators obtained in any such GIE are gauge invariant by construction, so there is 
actually a plethora of physically (and numerically) inequivalent gauge invariant decompo- 
sition schemes. 2 Still, the GIE based on the LC gauge ]l^, [l(], ||] has a quite distinguished 
status because it is the only GIE relevant to high energy experiments, as evidenced by the 
fact that AG obtained in this particular GIE (but not anything else) coincides with the 
experimentally measured gluon polarization. Therefore, when restricted to the context of 
high energy QCD, practically we do have a complete, gauge invariant decomposition of the 
nucleon spin. By introducing the so-called potential OAM 

L pot ^L"-Ll an , (1.4) 

one can make the following connection between (fb^) and (a GIE of) (|1.3|) 



Ji = ±AV + Ll an + L pot , (1.5) 



J° + L pot = AG + L c ' an . (1.6) 



In the case of the LC-gauge GIE, all the entries in ( |1.5| ) and (L6) can be explicitly written 
as the matrix element of manifestly gauge invariant operators Q. 

The goal of this paper is to understand ( |1.5| ) and ( |1.6| ) at the density level, L q can = 
j dxLcan(x), etc., where x is the momentum fraction of partons. We shall derive exact 
expressions for Ltan(x) and L 9 ca n{x) in terms of the twist-two and newly defined twist-three 
GPDs. Via the equation of motion, the latter can be written as the sum of the twist-two 
part ('Wandzura-Wilczek contribution') and the 'genuine twist-three', quark-gluon and 
three-gluon correlators. These results fit well with the standard framework of perturbative 
QCD, and therefore serve as the starting point for further first-principle calculations. We 
finally show the ^-moments of the results and confirm that the lowest moment reproduces 
(O) and Q. 



Contrary to the claims in |M [UJ, h2| . See |15| for the latest discussion. 
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2. Quark canonical OAM 

We start by defining the nonforward proton matrix element of the quark-gluon ('F-type') 
twist-three operator 3 

^^ e i|(*i+aa)+iAi(*a-*i)(p'^'|^(_An/2)7 + W_A xip(Xn/2)\PS) 

^TT 27T 2 2 

= \p + e +i »°u{P'S') lblp u(PS)A (T <S> F (x 1 ,x 2 ,^t) + ■■■ , 

^P + e +i ^S p A a $ F (x u x 2 ) + --- , (2.1) 

^^ e i7(»i+^)+to-=i)(P'5'|^(_An/2)757 + W_A gF +i (fin)W ^(Xn/2)\PS) 

= - l -P+e+^u(P'S') l5lp u(PS)A^ F (x l ,x 2 ,^t) + ••• , 

„ -iP+e +lpa S p A^ F ( Xl , x 2 ) + ■ ■ ■ , (2.2) 

where P' 4 ~ 5^P + is the proton momentum taken in the infinite momentum frame in the 3— 
direction, and S 1 ^ is the longitudinally polarized spin vector normalized as P 2 = —S 2 = m 2 , 
the nucleon mass squared. We have also introduced the average 

(P + ,0, 0,P~) and a lightlike vector = 5^/P + . The Wilson line W makes the nonlocal 
operators gauge invariant. 4 The momentum transfer is denoted as A^ = P' p — P p from 
which we define the skewness parameter £ = —A + /2P + and the Mandelstam variable 
t = A 2 . We assume that A^ 1 is small, and in the above expansion we have kept only the 
linear term in A which contains both the structure ~ e^Aj and a factor of S + = S 1 -. 5 The 



latter is apparent in the third line of ( |2.lD and fl2.2[ ) where we approximated ~ S ,fl in 
the spinors (admissible due to the explicit factor of A CT ) and denoted this common vector 
as S^. By the same token, we neglected the dependence on £, t in <& F and <& F (and in 
all the twist-three distributions defined later). In the following, we further assume that 
A + = 0, namely the proton is elastically scattered. This turns out to be a convenient choice 
because, when calculating the OAM, one differentiates nonforward matrix elements with 
respect to A %=1,2 (but not A + ) and then takes the limit A M — > (see, e.g., (^j|) below). 
Thus one can actually set A + = from the beginning. As we shall see, in the discussion 



3 Our conventions are 75 = — i7°7 1 7 2 7 3 , e 0123 = +1, and g^ v — diag(l, — 1, — 1, — 1). The light-cone 
coordinates are defined as x = -7= (a; ± a; 3 ) = x=p. The four-vector indices are denoted by Greek letters 
a, jit, ... and the transverse coordinates (momenta) are denoted by the indices i, j = 1, 2. We also introduce 
the two-dimensional antisymmetric tensor = — e lJ = — ejj , e 12 = +1. 

4 We shall use the same notation W both for the fundamental and the adjoint Wilson lines. The subscripts 
of W indicating the initial and final points will be often omitted when they are obvious from the context. 

5 We note that the spinor product in ( |2.l[ ) and fl2.2|) can be written in a seemingly different way using 
an identity 

ie +tp ' 7 u(P'S') l5 'y p u{PS)A„ = 2P + u(P'S')Yu(PS) , (2.3) 

which follows from the Dirac equation = {f — m)u(PS) = u(P'S')(i*' — m). The linear dependence on 
A (and S + ) is not manifest in this alternative form. 
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of the gluon OAM this choice is also very helpful when relating the matrix elements of 
different operators via the equation of motion. 

From PT-invariance, we find the following symmetry properties 

$ F (xi,x 2 ) = -<^ F (x2,x 1 ) , $ F (xi,x 2 ) = $f(x2, Xi) . (2.4) 

As shown in the function <&p(x\,X2) is related to the potential angular momentum 
(0) 

L pot = I dx 1 dx 2 V <£> F {xi,x 2 )= [ dXdx V- $>f(X, x) , (2.5) 

J Xl-X 2 J X 

where V denotes the principal value, and in the second equality we switched to the notation 
X = Zi±^, x = Xl -x 2 . 

Next define the l D-type' twist-three distributions 

2lT 2-7T 2^ ^2 

^e + ^S p A^ D (x 1 ,x 2 ) + -- - , (2.6) 
— ^e^+^+^-^HP' S'\^{-\n/2) lbl + W 

27T 2tt 2 2 
Kie +i ^S p A a $ D (x 1 ,x 2 ) + -- - , (2.7) 

where ^ a = \(D a - lJ a ), D a = d a + igA a , t> a = *d a - igA a . The structure of the 
right-hand-side has been fixed by the same criteria and approximations as in (|2,ip and 



( |2.2j ). This time we have 

§d(xi,X2) = $d(x2,xi), &d(xi,x 2 ) = -&d(x2,xi). (2.8) 
The integral of <£> D gives Ji's quark OAM L q |H§ 

L« = f d Xl dx 2 $ D (xi,x 2 ) = -ir lim (P' S'\^(Oh + i^^O)^(0)\PS) . (2.9) 

J 2b + A->0 io/l\ l 



There are generic relations between the F-type and D-type distributions |19j which 
can be derived from the following identities 



i 



W_x u D l (nn)W u x = — r dtW_x t gF +i (tn)W t x + W^xxD\Xn/2) 

* .1 A 2 

2 



W_x, SV)^, A = -ttt r dt^_ A , 5 F +i (tn)Ty, A + < D\-\n/2)W_xx . (2.10) 

2 ^ ^2 P + J X 2 2 22 

Using these identities, we find 

$d(x!,x 2 ) = V *f(jci, x 2 ) + S(x! - x 2 )L q can (xi) , (2.11) 

Xi — x 2 

$d(xi,x 2 ) = V $ F (x u x 2 ) , (2.12) 

X\ - x 2 
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where 




X ■ 



1 (W_x xiDHXn/2) - i5H-\n/2)W_x x) ^{\n/2)\PS) 

2 V 2 2 2 2 / 

" 2~h I (e(/u " A/2) + e(/i + A/2)) 

x(P'5'|^(-An/2)7 + VF_A^ 5 F + ^/in)W^AV'(An/2)|P < S}| . (2.13) 

(e(x) = x/\x\ is the sign function.) In ( |2.12| ), the delta function term 8{x\ — x?) is not 
present because $_d is antisymmetric in x\ and x-i- As already implied by the notation, 
the integral of L% a n gives quarks' canonical OAM 

L\ an ee J dxLl an (x) = ^ km -A_^(P'S^(0) 7 +f% re (0)^(0)|P ( S) , (2.14) 



where 



D% ure (Xn) = D a + X J dfiKQi - X)W x ,F +a ( fm ) 
= D"-igA« hys (\n) 

^d a + igA« ure (\n), (2.15) 

with /C(/i) being either |e(/i) or ±0(±^), 6 is the gauge covariant generalization of the 
canonical momentum id a in the LC-gauge GIE [fTPf . ( |2.15| ) also defines the 'physical' and 
'pure gauge' parts of the gauge field A a = A^ hys + Ap Ure which play a key role in the 
construction of || Q. 7 We can thus regard fl2.11| ) as the doubly-unintegrated version of 
the relation L q = L q can + L pot corresponding to the decomposition iD a = iDp Ure — gA^ hys . 
We see that at the density level, the decomposition of L q into the canonical and potential 
parts is very natural. 

Note also that, due to the delta function in fl2.11|) , the canonical OAM density has 
support only at x\ = x%, namely, the gluon has zero energy and the outgoing and re- 
turning quarks have the same light-cone momentum in the quark-gluon correlator. This 
implies that the argument x of L% an {x) can be indeed interpreted as the momentum 
fraction of quarks, in much the same way as in the usual parton distribution function 
(PDF), and also suggests the uniqueness (or the preferred choice) in defining the density 
Lean = J dxLcan{x). On the other hand, from this point of view, there seems to be an 



6 Note that, in Eq. (2.13), one can replace |e(^ + A/2) = ±(6(±(ji + A/2)) - ±) with ±0(±(fi + A/2)) 
(similarly for |e(/i — A/2)) because the /i-integral for the '±s' term is unrestricted, and this corresponds 



to the point X\ = X2 in (2.1) where the function <1>f vanishes. The different choices of K. correspond to 
different boundary conditions in the light-cone gauge. See jTTj] for the detail. 

7 The authors of |j, [?| work in the GIE based on the Coulomb gauge, and accordingly their A p h ys is 
different from ours. 
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ambiguity in denning the density of the dynamical OAM (|2.9| ), L q = f dxL q (x). 8 [See, 
also, the discussion in [18].] 

Let us derive a relation between the OAM and the generalized parton distributions 
(GPD). From the equation of motion, one can show that 

tP{-z-/2)f(WD_ - D ^W)tp(z- /2) = i)^ + {WD l - frwty 

+ie ij ^ 5 -fj{WD^ + 15-W)ip - ie ij ^ 5 -/ + {WDj + ^5jW)tp . (2.16) 

(In our convention, L>_ = D + .) The nonforward matrix element of ( |2,16| ) reads 



2- 







(P'S'\i>{-z/2)YW^(z/2)\PS) = {P'S'\{jj + (WD i - frwtylPS) 



dz 

+2ie ij V^{P'S'\i)^jW^\PS) - ie ij {P'S'\i)^ + {WDj + DjW)^\PS) .(2.17) 

The second term on the right-hand-side contains the total derivative T> corresponding 
to the translation of spatial coordinates: V^O{— z, z) = lim a M_>.o z + a, z + a) — 

0(—z,z)). This term vanishes since, in the nonforward matrix element, T>^ ~ i[P'^ — 
P^) = iA^ and we assume A + = 0. The matrix element of the left-hand-side may be 
parameterized by the twist-two [jlj and twist-three M GPDs 

{P'S'm-z/2)^Wi;(z/2)\PS) = J dxe- ixP+z ~ {(#,(*,£,*) + E q (x,£,t))u(P'S')^u(PS) 

-E q (x^,t)—u(P'S')u(PS) + Gsix^^tMP'S'^uiPS) + ■■■ ,(2.18) 
m J 

where 7^ = S^j 1 . In (p.180 , we have kept only the terms relevant to us. 9 For the same 
reason as explained below ( |2.2| ), the arguments £ and t of GPDs will be neglected in the 
following. Noting that (cf. (pip 

u(P'S')fu(PS)^ie^A :j ^ T , (2.19) 
and comparing the terms proportional to e lJ AjS + , we obtain 

x{H q [x) + E q (x) + G 3 (x)) = H q (x) + ^Jdx' ($ D (x, x') + $ D (x', x) + <b D (x, x') - $d(x', x) 

= H q (x) + L\ an (x) + J dx'V^-^ ($ F (x, x') + $ F {x, x')) , (2.20) 

where H q is another GPD which gives the polarized quark distribution H q (x,A = 0) = 
Aq(x) in the forward limit. Integrating over x, we find, using (2^) and ([2.5|), 

Aq + L q can + L pot = j x(H q (x) + E q {x) + G 3 (x))dx . (2.21) 

8 For instance one may define L q (x2) = J dx\<br>{x\, X2), or L' q (X) = J d{x\ — 22) $d{x\, X2) giving 
different functions L q (x) ^ L' q (x) in general. 

9 There are different conventions for the parametrization of twist-three GPDs in the literature (cf., 
Refs. ^0], ^|). Here we take any one of such parametrizations, extract the structure u'-f^u (by using 
(p.3|), if necessary [H]), and redefine its coefficient to be 6*3(0;). 
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We now recall that J q = \ f x{H q (x) + E q (x))dx = L q + \t\q. Then it follows that 

J dxxG 3 (x) = —L q . (2.22) 

This formula was previously derived in albeit in an approximation ('parton model') 
which neglects gluons. 

To proceed further, we need another independent relation which can be used to elimi- 
nate G3 from ( [2.20 ). For this purpose we employ the identity pi] 



—^(-z/2h a WiP(z/2) - —$(- z /2) lfX WiP(z/2) 



2 



■z^ p V» {^(-z/2) l5 YWi;(z/2)) 



/ \ due a ^{-z/2) lb yWgF^{uz)z T WiP(z/2) 



iz^ / duuii;(-z/2)-f fl WgF aT (uz)z T W^(z/2) . (2.23) 



In the above nonlocal operators, z 11 is generic, not necessarily proportional to n^. The 
nonforward matrix element of the right-hand-side reads, in the case of z^ = 5^Lz~ and for 
the component a = i, 

^-e^A,(P'S'|Vi(-z/2) 75 7 + ^(z/2)|P5) 



(z-) 2 12 



2 /_i 

2 
1 



dr (P'S'\4>{-z/2) 1 ^ + WgF +l {rz)W^{z/2)\PS) 



+i{z-) 2 r dTT{P'S'\$(-z/2)-/ + WgF +l (Tz)Wip(z/2)\PS) . (2.24) 

"2 

The first term is related to the polarized quark distribution 

^A 3 (P'S'\^(-z/2)^ + W_ u i;(z/2)\PS) = -^A,£^ J dx e~ lxp+ z ~ ^H q {x) , 

(2.25) 



while the second and the third terms become 

(z-) 2 n 



dT{P'S'\^(-z/2)j 5 j + WgF +t (Tz)W^(z/2)\PS) 



(2.26) 



i( z T J x dTT{P'S'\i>{-z/2) 1 ~ ] -WgF +l {Tz)W^{z/2)\PS) 

= -ie"A^ J dxie -^~ J dX2 v^- (A - A) », (XI ,„). ( 2 . 27) 
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Next consider the matrix element of the left-hand-side of ( 2,23; ). Away from the 
light-cone, ( 2.18[ ) may be covariantly generalized to 10 

(P' ' S'\^{-z/2)^Wij{z/2)\PS) 

= Jdx e- ixp z u(P'S') ({H q + E q + G 3 ) 7 M - G 3 ^-P^ u(PS) + ■ • {2.28) 



This leads to 



JL{P>S'\iP(-z/2)fWiP(z/2)\PS) - ^(P'S'\^(-z/2)^WMz/2)\PS) 

J dxe -™P + *- (±( x (H q (x)+E q (x)))+x-^G 3 (x)\ , (2.29) 

where we have set z^ = S^z" after differentiation and kept only the 0(A) terms using 
(|2~19|) . We thus find 

±(x(H q (x) + E q (x)))+x^G 3 (x) = ^H q {x)+ I dx'V-^ ( £ - $ F (x, x') 
ax 'ax ax J x — x \ ox ox J 

( 2.30 ) can be formally solved for G3 with the boundary condition G 3 (x = ±1) = 0. Sub- 
stituting the result into ( 2.20| ), and performing integration by parts, we arrive at 

Llanix) = x / ^r(H q (x') + E q (x')) - x / ^H q (x>) 

Jx X J x X 



f-e(x) rl 

X / dx\ / dX2&F(xi, X2)V- 
Jx J-l * '- 



2>X\ — X2 



x\{x\ - X2) 2 

r-e(x) rl ^ 

-x / dx\ / dx 2 $F(xi, X2)V—~-, r. ( 2 -31) 

Jx J-l Xf{xi-X 2 ) 

This is an identity which relates the quark canonical OAM with the twist-two GPD and 
genuine twist-three quark-gluon correlators. By substituting ( |2.31[ ) into ( |2.11| ), one finds a 
similar identity for the density <&d(xi,X2) of Ji's quark OAM (cf., (|2.9| ), see, also footnote 
8). In the parton model ('Wandzura-Wilczek approximation'), one neglects the $-terms 
after which there is no distinction between Lean and L g . 

Finally in this section, consider the moments of ( |2.31 ). Using the formula 

f ±:L f ±1 dz 1 r ±l 

J dxx n ~ l J —f(z) = -J dxx^fix), (2.32) 

we find 

/ dxx n - l L q can {x) = — *— / dxx n (H q (x) + E q (x)) l — f dxx n - l H q (x) 

J-i n + 1 J-i n + 1 7-i 

l — C dxdx'x n - x (p 3x - f q> F ( x , x >) +V—!—$ F ( X:X i)) . (2.33) 

n + 1 7-i \ {x — x ■ y x — x J 



10 Note that the GPD E q should be kept in the fj, = + component because we are interested in the 0(A) 
terms. 
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In particular, when n = 1, 

Llan = j dxL q can {x) = \J dxx{H q (x) + E q (x))-^ J dxH q (x) 

If / 2>x — x' 1 \ 

-- / dxdx'iV-. -2$ F (x,x') + V -$ F (x,x')\ 

2 J \ (x — x'Y x — x' J 

= J q - ^Aq - L pot , (2.34) 



where we used ( [2.5| ) and the symmetry properties fl2.4|) . ( [2.34 ) is the correct spin decom- 
position formula (|1.5|) in the quark sector. 



3. Gluon canonical OAM 

The analysis of the gluon canonical OAM is entirely analogous to the quark case. We first 
need several definitions. 

• Unpolarized/polarized gluon distribution functions: 
^ e iXx (PS \F +a (0)WF + P (An) I PS) 
= -\xG{x)(P + f(g aP - P a nP - P?n a ) - % -x/\G{x)P + e + ~ a ? S + + ■■■ . (3.1) 
Gluon GPDs: 



d\ 1 / \ n a " 

^e iXx {P'S'\ - - ( F aT {-Xn/2)WF p T {Xn/2) + F pT WF a r ) + ^—F pT WF pT \PS) 
2tt 2 V / 4 

1_1 p( a 'jrr^ T A 
= -H g (x)u(P'S')P^u(PS) + -E g (x)u(P'S') — ^u(PS) 

+ l -F g {x)u{P'S')P^u{ p S) + • • • , (3-2) 

where A^B^ = ^(A a B 13 + A^B a ). In addition to the usual twist-two GPDs H g 
and E g , we have introduced a twist-three GPD F g which satisfies J dxF g (x) = 0. 11 
There are other twist-three terms, but those are not needed for the present purpose. 

F-type three-gluon correlators: 

_1 r ^d^ e ih( xl+X2 ) +ifl ( X2 - xl ), pl S '\F +i {-Xn/2)WgF +j (fin)WF +k (Xn/2)\PS) 
[P+Y J 2n Z7r 

« (^M(x 1 ,x 2 )g ik e + ^ + M(x 2 ,x 2 - x 1 )g i 'e+ k ^ - M(x u xi - x 2 )g jk e +i ^)s p A a + ■■■ . 

(3.3) 

In order to fix the tensorial structure of the right-hand-side, we followed the method 
developed for the three-gluon correlator in the transversely polarized case p2|, 23]. 



11 After the ^-integration, the operator on the left-hand-side of ( |3.2| ) becomes a local operator. Its matrix 
element should not depend on non-covariant expressions such as and Vf- 
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[As in (|2.lD , S p should more precisely read \u{P' ' S')^^ p u(PS).} By PT invariance, 
it follows that 

M(x u x 2 ) = -M(x 2 ,x 1 ). (3.4) 
On the other hand, from permutation symmetry we find 

M(xi,x 2 ) = M(-x 1 ,-x 2 ). (3.5) 



Contraction of i with k leads to an expression similar to (2.1) 

f ^^ e ^ (a;i+X2)+ ^ (x2 - xl) (P , 5 , |F +Q (-An/2)iy 5 F +J '(^n)VFF + Q (An/2)|P5) 

J 2tT 2tT 

« (P + ) 2 e +jpa S p A a M F ( Xl ,x 2 ) + ■■■ , (3.6) 
where we defined 

M F {xi,x 2 ) = 2M(x 1 ,x 2 ) + M(x 2 ,x 2 - x x ) - M{x 1 ,x 1 - x 2 ) , (3.7) 

with the property Mp(xi,x 2 ) = —Mp(x 2 ,xi). We also note the expression in the 
coordinate (z~) space 

(P' S'\F +a ((n)WgF +j (fin)WF + a (\n)\PS) 

« {P + ) 2 e +jpa S p A a J dxxdxie-^-^-^+^Mpix^x-z) + ■■■ . (3.8) 

D-type three-gluon correlators: 

^^e i T(*i+^)^(*a-*i)/*(P'^qF +a (-An/2)W r ^ i (Mn)Wi^ .(An/2)|P5) 
« P+ e +^5 p A CT M Z) (x 1 , x 2 ) + • • • , (3.9) 



with Md(xi,X2) = M£)(a;2,xi). One may define, as in ( |3.3|) , more general D-type 
correlators in which the Lorentz index a is not contracted. However, we do not need 
them for the present purpose. 

In order to make contact with the gluon OAM, we need three-gluon correlators in 
which one of the F's is replaced by A p hy S . They can be easily obtained via the convolution 



(P'S'\F +a {(n)WgF +l (fin)WA p a hys (\n)\PS) 

« P + e +ip,7 S p A a J d Xl dx 2 e~" lA -^ 2 - Xl ^ +iX2C i/C(xi)MHxi, x 2 ) + • • • ,(3.10) 



{^P>S'\F +a {i:n)W^\^n)WAf ys {Xn)\PS) 



e+ ip,J S p A a J dxidx 2 e- ixiA - i(x2 - Xl) ^ 2 ^/C(xi)M D (xi, x 2 ) + ■■■ , (3.11) 
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where K{x) = V- in the case /C(/z) = ^e(/u) and K.(x) = in the cases K.(/j,) = ±0(±/z). 
As in ( |2.1lD and (|2.12| ), the F-type and D-type distributions are related: 

M D {xi,x 2 ) M F (x 1 ,x 2 ) xf , Tq , , ,„ 10 , 
= — -5 (xi -x 2 )L 9 can (xi). (3.12) 

X\ X\{X\ - x 2 ) 

The function L 9 can (x) thus defined can be identified with the density of the gluon canonical 
OAM f dxL 9 can (x) = Lean because it satisfies [Q] 12 



ie + ^S p A a j\xL 9 can (x) = -(P'S'\F +a ^; ure AP h y s \PS) . 



(3.13) 



Moreover, as already argued in the quark case, the delta function 8{x% — x 2 ) suggests the 
uniqueness (or the preferred choice) of the density L 9 can = f dxL 9 an (x) with the variable x 
interpretable as the momentum fraction of gluons. 



Let us find the counterpart of ( [2, 16 ) in the gluon case. Consider the following nonfor- 
ward matrix element 

/ = -^-{P'S'\F +T (-z- /2)WF i Jz-/2) + F iT (-z- I2)WF + Jz- I2)\PS) . (3.14) 

oz~ 

Since A + = by assumption, we can freely make translation of the field coordinates in the 
z~ direction and obtain 

/ = (P'S'\F +T {0)WD^F i T (z-) - F iT (-z-)D_WF + T (0)\PS) 
= (P'S'\F +T (0)W (piF+^z-) - p) T F + \z-)} 

+ f-F+ T (-z-0 + F +i (-z-)%) WF +T (0)\PS) , (3.15) 



where in the second equality we used the Jacobi identity. We then use the formula 

-F +T (0)W~3 T F +i (z-) = F +T l5 T WF +i - Vj (F +j WF +l ) 

-i [ duj- F +j {0)WgF + Auj-)WF +i (z~), 
Jo 

and the equation of motion D T F^ + = gip'y+^ip to get 



IP'S' 



F +T (w~B' - frw) F + T + Vj (-F + iWF +t + F +i WF + i) (3.16) 



-il)(-z-)^ + WgF +l {0)W4>{-z-) + ^(z-)-y + WgF +i (0)W4>(z~) 

r*~/ 2 / . .\ 

+i J cLj- [-F +1 {-z- /2)WgF + j (uj-)WF +l (z- /2) + F +t WgF+ j WF +3 ) 



\PS) 



12 There is an overall sign error in the definition of the gluon canonical OAM in Q|. The one-sided 
derivative Dp Ure in the definition of L c g an can be replaced by the symmetrized one Dp Ure without affecting 
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These matrix elements can be expressed by the three-gluon correlators and AG(x). On 
the other hand, the starting expression ( 3.14| ) can be directly expressed by the GPDs. By 



equating the coefficient of S + e l ^Aj in these two expressions and using ( 3.12| ), we find 



\{H g {x) + E g (x) + F B (x)) - AG(x) + 2 J dX * F &' x) - 2L° an (x) 

f , ,^M F (x,x') f , .M F (x,x') 
= -2 / dx'V P ' - 2 / dx'V ; v ' 

J x{x — x ) J x{x — x ) 

= -4 / dx> v M(x,x>) + M(xW- X ) (317) 
J x(x — x') 

where we defined a function 

M F (x, x') = M(x, x-x') + M(x', x' - x) . (3.18) 

Note that Mp(x,x') is symmetric under x x' , although M(x,x') and M F (x,x') are 
antisymmetric. We now integrate ( p.17 ) over x, recalling the basic properties 



~J dx(H g (x) + E g (x)) = 2J g , J dx AG '(z) = 2 AG , J dxF g (x) = , (3.19) 



as well as (|2.5|). The result is 



n f, r ,M(x,x') + M(x',x' - x) 

J g - AG + L pot - L 9 can = -2 / dxdx' V y ' ' . K -£ - . 3.20 

J x(x — x') 



Using the antisymmetric property ( |3,4| ) of M(x,x'), one can easily check that the right- 
hand-side of ( |3.20| ) is identically zero. Eq. ( |3.20[ ) is then nothing but the spin decomposition 
formula ([Eq). 

Next we consider the counterpart of ( |2.23| ) in the gluon case. The generalization of 
fly) off the light-cone z» ^ 5tz~ is (cf.,(g2§) 



;^<P / S / |P QT WP^ r + P^P a r |PS') = - /" dxe- ixp - z u{P'S') 

x ((P • z 7 a + P a z • 7)(fl" ff (x) + P 9 (x) + F g {x)) - 2P a 7 • zF g )u{PS) + ■ ■ -(3.21) 

This immediately gives 
d 



(P'S'\F iT WF p T + F /3T WF i T \PS) - (i m) 



« £L| dxe-^" |A( x(jffg ( x)+jBa(a , )) ) +x 2^_EM\ ie v Aj S+, (3.22) 

where again we set z^ = 5^_z~ after differentiation and kept only the linear terms in A. On 
the other hand, the left -hand-side of ( |3.22j ) may be computed directly using the equation 



- 12 - 



of motion and the Jacobi identity 

^ {jz~^ z P( P ' S '\ FiTWFP T + F Pr WF i T \PS) -(if) n)\ 

= {z-f {P' S'\\ -gi){^ + W F +i {z)W^{Q) + <7^(0)7 + W F +i (-z)W^(<d) 



-Vj (F +j WF +i ) + Vj (F +i WF +j ) 

/•1/2 

+iz~ / du (F +l WgF + Auz)WF +1 - F +3 W gF + Auz)WF +l 

J -1/2 V 

1/2 



+2iz- j duuF +T WgF +i (uz)WF + T ^\PS) . (3.23) 
[Again, = 5^_z~ after differentiation.] The matrix elements in ( [3.23 ) can be evaluated 



as 



z~ [ dxe~ ixp+z ~ i-2 [ dX-^-$ F (X,x) + — (xAG(x)) 
J I J ox dx 

Equating this with ( |3.22| ), we obtain 

||(» WW + * s W)) + ^3£> = _ 2 / d x J-* F <*, x) + A (lAG(l)) 

- 2 /^rb(s + il)*'^ 



which can be solved for F g (x). Eliminating F g from ( 3.17 ) in this way, we finally arrive at 

r M 35 ) rfr' rfr' 
^n(^) = ? / 4(^,00 + ^0) " x / 32 AG(x') 

Z Jj. X J x X 

re{x) faj r re(x) rl ]_ 

+2x J — J dX$ F (X,x') + 2x J dxi J dx 2 M F (x h x 2 )V 



x\{x\ — X 2 ) 

fff 1 ) r 1 2xi — xo 
+2x d Xl dx 2 M F {x 1 ,x 2 yP- sr ^ (3.26) 

Jx J — 1 X 2 ) 

This is our main result which expresses the gluon canonical OAM density in terms of 
twist-two GPDs and genuine twist-three contributions & F , M F and M F . One may wish 
to neglect the latter in the spirit of the Wandzura-Wilczek approximation. In that case it 
may make more sense to neglect only the last two terms M F and M F since they do not 
contribute to the integrated gluon OAM Lean = j L 9 can {x), see below. 
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In the moment space, ( |3,26 ) takes the form 



J dxx^I^ix) = I dxx n -\H a {x)+E a {x)) - -L^J dxx n - l AG(x) 

+ ^Tl J dxxn ~ 2 J dX$ F (X,x) 

+ f dxdx'x n - 2 V MFix,X,) +— ( dxdx'x n - 2 V 2X ~ X ' M F (x, x') . (3.27) 

n + 1 J x — x n + 1 J (x — x ) z 

In particular, when n = 1, we get 

LLn = J 9 + L pot -AG + J dxdx'V ^^ + J dxdx'V x 2x 2*' l)2 M F (x,x') 

= J g + L pot — AG + f dxdx' v M F ^ X ')+M F (x,xi) (3 28) 

J x(x X ) 

where we used the antisymmetric property of Mp. As already noted, the last integral 
vanishes and we recover the sum rule (Of)- 

4. Conclusions 

Due to the difficulty in experimental measurements, twist-three GPDs have received rela- 
tively little attention so far, and even less attention in the particular context of longitudinal 
polarization. However, we have seen that they are intimately related to the decomposition 
of the nucleon spin beyond (|l.l| ) , and therefore provide strong motivation for further study. 
Especially it is an urgent task to pin down experimental processes which are sensitive to 



these distributions. In this respect, the connection [24, [4| between the canonical OAM 
and the OAM from the Wigner distribution (constructed from the 'generalized parton 
correlation function' p5[) may be helpful. 
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